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We have performed measurements of scintillation light in hquid neon, observing a signal yield 
in our detector as high as (3.5 ± 0.4) photoelectrons/keV. We measure pulse shape discrimination 
efficiency between electronic and nuclear recoils in liquid neon from 50 and 300 keV nuclear recoil 
energy. We also measure the CaS parameter in liquid neon between 30 and 370 keV nuclear recoil 
energy, observing an average CcS= 0.24 above 50 keV. We observe a dependence of the scintillation 
time distribution and signal yield on the pressure and temperature of the liquid neon. 

PACS numbers: 61.25.Bi,25.40.Dn,34.50.Gb 



I. INTRODUCTION 

The DEAP / CLEAN collaboration has proposed to use 
liquid argon and hquid neon as scintillating targets in 
large scale detectors to observe both pp-solar neutrinos 
and dark matter in the form of weakly interacting mas- 
sive particles (WIMPs) [IHli- In particular, the Mini- 
CLEAN experiment is designed to do a dark matter 
search with first liquid argon and then liquid neon as both 
a systematic check on backgrounds and as a precursor to 
a multi-tonne scale detector 4] . Because of the expected 
v4^-dependence of the cross section for WIMP-nucleon 
scattering (5|, filling the same detector alternately with 
argon and neon will change the expected dark matter sig- 
nal in a known way while maintaining an identical level of 
external backgrounds. Here, we extend previous studies 
of the scintillation properties of liquid neon Q . 

Like the other liquefied noble gases, liquid neon is rel- 
atively inexpensive, is easily purified of radioactive con- 
taminants, scintillates brightly when exposed to ionizing 
radiation and is dense enough to self-shield, reducing the 
background level in the center of a larger volume of liq- 
uid. The key for dark matter detection is to be able to 
suppress the electronic recoils that make up most of the 
backgrounds from the nuclear recoils that would make 
up a WIMP signal by use of some combination of self- 
shielding and discrimination. Because Ve — ^ scattering 
also produces electronic recoils, discrimination does not 
improve the performance of a neutrino detector. How- 
ever, neon has no long-lived radioactive isotopes and is 
more readily purified of such contaminants as '^^Ar and 
^^Kr 0] than the heavier noble gases, rendering it an 
ideal target for low-energy neutrinos. 

When ionizing radiation interacts in liquid neon or any 
other liquefied noble gas, ultraviolet light is produced 
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via scintillation. The incoming radiation collides with an 
electron or nucleus in the liquid and deposits energy. The 
resulting electronic or nuclear recoil then excites or ion- 
izes other neon atoms. The excited atoms quickly com- 
bine with surrounding ground state atoms to form dimer 
states, Ne2. The ions also bond with ground state atoms 
to form ionized molecules, NeJ, which can in turn recom- 
bine with free electrons to also form the excited dimer 
states. Finally, the dimer states decay emitting scintil- 
lation photons with wavelength 80 nm. Because these 
photons have insufficient energy to excite neon atoms [1] , 
the liquid does not absorb its own scintillation light and 
large light yields can be attained from liquid noble gas 
detectors. 

The metastable molecules are produced in both singlet 
and triplet states. The singlet state decays in nanosec- 
onds, while the triplet state undergoes a forbidden spin 
flip before decaying, extending the triplet lifetime to 
15 /^s in neon [6|. Because nuclear and electronic re- 
coils produce different ratios of singlets to triplets, pulse 
timing can be used to identify the source of the initial 
excitation via pulse shape discrimination (PSD) [2], 
As the PSD depends directly on the scintillation timing, 
we have attempted to re-measure the time constants of 
liquid neon scintillation in this work. We discovered two 
additional time components in the light signals produced 
by ionizing radiation, and we also found that both the 
intensity and lifetime of each component changes signif- 
icantly with the temperature and pressure of the liquid, 
a new feature that is not observed in the heavier noble 
gases. 

A second property of scintillation in liquefied noble 
gases is that nuclear recoils produce less light than elec- 
tronic recoils of the same energy. The ratio of signal 
yields at zero electric field for the two event classes is 
known as the nuclear recoil scintillation efficiency, or £off- 
The value of £cff in combination with the discrimination 
efficiency sets the nuclear recoil energy analysis thresh- 
old of a liquefied noble gas dark matter detector. Be- 
cause in general £cff < 1, two energy scales may be em- 
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ployed, denoted "keVee" and "keVr." The actual energy 
deposited by a nuclear recoil is expressed in units of keV 
recoil or keVr. Most calibration sources used in liquid 
noble gas detectors are 7 or /? sources that produce elec- 
tronic recoils in the liquid. Therefore, the energy cali- 
bration gleaned from these sources refers to the energy 
deposited by an electronic recoil, or keV electron equiva- 
lent (keVee). In practice, we often have no direct energy 
calibration for nuclear recoils and so use electron equiv- 
alent energies. The conversion factor is C^s- At zero 
electric field 



£;[keVee] = E\keVv\ x 



(1) 



PSD has been studied in argon [o'-TT'l and xenon [T^ 
Our group performed a measurement of PSD in neon 
as welL but nuclear recoils were observed at only one en- 
rgy 0- Similarly, C^s has been measured in argon 15 1 



and xenon [16li2(J|. but at only one energy in neon. In 
this paper, we report measurements of scintillation light 
in liquid neon over a wide energy range. We measure 
PSD from 50 to 300 keVr and C^s from 30 to 370 keVr . 
We also present results on the time dependence of liquid 
neon scintillation light as a function of temperature. 



II. EXPERIMENTAL APPARATUS 

The apparatus used in these measurements is named 
MicroCLEAN and has been described in detail else- 
where jSi]. A schematic is shown in Fig. [TJ The de- 
tector has an active volume of 3.14 liters viewed by 
two 200-mm-diameter Hamamatsu R5912-02MOD pho- 
tomultiplier tubes (PMTs), which are specifically de- 
signed for use in cryogenic liquids. The active volume 
is defined by a polytetrafluoroethylene (PTFE) cylin- 
der 200 mm in diameter and 100 mm in height, with 
two 3-mni-thick fused-silica windows at top and bot- 
tom. All inner surfaces of the PTFE and windows are 
coated with (0.20 ±0.01) mg/cm^ of tetraphenyl butadi- 
ene (TPB) f2l|, which shifts the wavelength of the ultra- 
violet light to approximately 440 nm. The TPB thick- 
ness on the PTFE cylinder was reduced by 33% rela- 
tive to that used in Q to mitigate TPB absorption of 
the wavelength-shifted light. The MicroCLEAN vessel 
is contained within a stainless steel vacuum dewar for 
thermal insulation purposes. 

The gas handling system is shown schematically in 
Fig. [5J Ultra-high-purity neon flows through a Nupure 
Omni 1000 22] 65j getter and an external charcoal trap 
at 77 K (not shown) before passing into a vacuum cryo- 
stat through a tube at the top. Previous studies have 
shown that activated charcoal provides a very effective 
purifier of many impurities in neon, including radon Q. 
Once inside the vacuum cryostat, the gas flows into a 
copper liquefler cell attached to the second stage of a 
Cryomech Model PT805 cryorefrigerator before dripping 
into the MicroCLEAN vessel. 
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FIG. 1: (Color online) Schematic representation of the Mi- 
croCLEAN scintillation cell. 



There is a circulation loop inside the vacuum dewar for 
continuous puriflcation of the neon, again using activated 
charcoal as the purifler. Liquid neon flows out of the bot- 
tom of the MicroCLEAN vessel into a nearby plumbing 
volume containing a heater. The heater (labeled "Circ. 
pump" in Fig. [2]) acts as a pump by boiling the liquid; 
the resulting pressure diflFerential between the heater and 
the liquefler drives gas up a tube through a charcoal trap 
mounted inside the cryostat before reentering the top of 
the liquefler, where gravity pulls the liquid back into the 
main chamber. 
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FIG. 2: The gas handling system for the MicroCLEAN runs 
described in this paper. The circulation pump consists of a 
heater that evaporates the liquid; the gas then flows through 
the charcoal trap before being reliquefied. 



The "83Kj-m Generator" shown in Fig. [2] refers to the 
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setup for introducing ^^Kt'^ into the liquid neon for cal- 
ibration purposes as described in Sec. Ill Bl and [23j . The 
83j^j.m generator is connected to the gas inlet line just 
outside the vacuum dewar. The incoming gas can be di- 
verted through the trap on its way into the detector to 
introduce ssj^j.m ^toms into the neon. 



A. Data acquisition and processing 

The data acquisition (DAQ) system consists of a 250 
MHz, 12-bit CAEN V1720 waveform digitizer (WFD). 
The two PMT channels are recorded with a record length 
of 64 /is to collect as much of the long-lived triplet light as 
possible, and a sample trace is shown in Fig.[3l The PMT 
signals are passed through a LeCroy Model 612AM dual- 
output amplifier, with one copy passed to the WFD and 
the other copy sent to triggering electronics. The trigger 
is set in several different configurations. In most cases, 
a trigger is generated when each PMT signal crosses a 
threshold of 1/2 times the mean height of a single photo- 
electron within a 100 ns coincidence window. For some 
low energy runs when an external detector is used as a 
coincidence tag, the trigger is generated when the sum of 
the two PMT signals crosses a threshold of roughly 1/2 
the size of a photoelectron in coincidence with a signal in 
the external tagging detector. In each 64 /is data record, 
approximately 8 /iS of presamples are collected to mea- 
sure the baseline of each pulse. The data collected by the 
DAQ software are saved in a ROOT-based file structure 
for analysis [24l|. 
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FIG. 3: Sample event from a single PMT digitized by the 
CAEN VI 720 waveform digitizer. More than eight microsec- 
onds of presamples are recorded to measure the baseline. 



performing this integral, we restrict the range of integra- 
tion to 100 ns windows around locations where the trace 
voltage crosses a threshold of three times the baseline 
root-mean-square. This integration method is a hybrid of 
single photoelectron counting and pulse integration. We 
reject data for which the trigger times of the two PMTs 
differ by greater than 20 ns. We define an asymmetry 
parameter. 



Sf 



St + -Sb ' 



(2) 



where St and Sb are the signal areas in the top and 
bottom PMTs. We reject data for which the asymme- 
try parameter is farther than 2a from the mean value as 
determined by fitting a Gaussian function to the asym- 
metry data. 

As the triplet component of liquid neon scintillation 
is spread out over tens of microseconds, there are many 
single photoelectron pulses well separated in time in the 
tails of events. To obtain a single photoelectron spec- 
trum for each PMT, the last 30 fj,s of each trace is evenly 
divided into 150 ns windows. If the voltage spread in 
the first and last 30 ns of each 150 ns window is less than 
three times the size of the baseline root-mean-square, the 
peak height and pulse area of the window are collected 
in a 2-D histogram; this restriction prevents inclusion of 
partial single photoelectrons in the data set. A Gaussian 
is fit to the 1-D histogram of the pulse areas to determine 
the mean single photoelectron size. 

In previous tests of the R5912-02MOD, it was found 
that the gain dropped substantially at 30 K [2^ . For the 
data presented here, we observe a sharp drop in PMT 
gain upon cooling to liquid neon temperatures, followed 
by a slow decrease in gain over several weeks before finally 
reaching a plateau. The low resulting gain necessitates 
the application of higher voltages to the PMTs and am- 
plification from the LeCroy 610AM amplifier; for most of 
the data presented here, the PMTs are biased at 1750 V, 
the absolute gain of the PMTs is w 1.1 x 10^, and the 
PMT signals are further amplified by a factor of 10 with 
the LeCroy amplifier. A single photoelectron peak is re- 
solved in the bottom PMT but not in the top PMT. To 
improve PMT performance, we use two high voltage lines 
to bias the PMTs, with a separate high voltage running 
directly to the first dynode. By tuning the first dynode 
voltage relative to the total voltage, we can achieve some 
improvement in the single photoelectron spectra but not 
enough to recover a resolvable single photoelectron peak 
in the top PMT. 



The data processing is similar to that described in [9j . 
The presamples are averaged to obtain a baseline and a 
baseline root-mean-square, and the baseline is subtracted 
from each trace. We integrate each trace from 50 ns be- 
fore the trigger (defined as the time at which the voltage 
first rises above 20% of its maximum value) to 44 fj,s after 
the trigger to determine the total area of each pulse. In 



B. Radioactive sources 

We use several radioactive sources for calibration pur- 
poses: the ®3j^].m generator described in [2^, a 10-/iCi 
^^Na source for 511-keV and 1275-keV gamma rays, 
and a D-D neutron generator producing a monoenergetic 
beam of 2.8-MeV neutrons [265. The ^^Kr"^ generator 
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produces metastable ^^Ki'^ nuclei that decay via two 
electromagnetic transitions with energy 32.1 keV and 9.4 
keV. The second transitions has a half life of 154 ns, 
and given the long time scale of liquid neon scintillation, 
both transitions are observed together as a single wave- 
form of total energy 41.5 keV. The ^^Na source produces 
positrons that immediately annihilate into two 511-keV 
gamma rays with opposite momenta; we tag the second 
gamma ray using a Nal crystal placed back-to-back with 
our apparatus to improve the data quality from the ^^Na 
source and reduce backgrounds. We reject data where 
the trigger times in the Nal crystal and the liquid neon 
differ by more than 30 ns, and we also make a cut to 
select events within 2a of the 511-keV peak of the Nal 
crystal. 

We also use a tagging detector when investigating nu- 
clear recoils with the neutron generator. We require a 
PMT viewing BC501A organic scintillator to record an 
event within 200 ns of an event in the liquid neon. The 
organic scintillator is contained in a cylinder with a diam- 
eter of 127 mm and a depth of 127 mm. We also make a 
PSD cut in the organic scintillator and a cut on the timc- 
of-flight (TOF) between the event in the liquid neon and 
the tagging detector. These cuts will be discussed in 
more detail in Sec. IIVI Both the neutron generator and 
the organic scintillator are located approximately 1.63 m 
from MicroCLEAN, and the experimental setup can be 
seen schematically in Fig. 01 By changing the scattering 
angle 6, we can choose the energy of the nuclear recoils 
that scatter just once in the liquid, Erec, using simple 
kinematics: 

-cos(6i)v/M2 + cos2(0) - 1], (3) 

where Ein is the incident neutron energy and M is the 
atomic mass of the target. We measure the angle by 
running a string marked at both ends and at its midpoint 
from the generator around the dewar and back to the 
generator. We mark the location of the midpoint of the 
string on the dewar and then repeat the process for the 
organic scintillator. The ratio of the arc length between 
the two points on the dewar to the circumference (175.6 
cm) provides the angle, 9. The accuracy is limited by our 
knowledge of the midpoint of the string, and by repeated 
measurements we estimate the uncertainty on the arc 
length to be 6.4 mm, for an uncertainty on the angle of 
1.3 degrees. 

C. Energy calibration 

We obtain an energy calibration by comparing ^^Na 
data with a simulation of the detector performed with the 
Reactor Analysis Tool (RAT) f27|. Combining Geant4 
and ROOT [3, [2^ in a single simulation package, the 
RAT simulation of MicroCLEAN contains a complete 
optical model of the inner detector and PMTs and has 
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FIG. 4: (Color online) Schematic of the neutron scattering 
setup. 



been used previously in studies of liquid argon scintilla- 
tion [l,[l^. In general, there are two features in the ^^Na 
data that can be used for calibration purposes: a 511-keV 
photoabsorption "hump" and the Compton edge. 

The analysis begins by first obtaining a rough estimate 
for the signal yield from the 511 keV hump, without using 
the simulation. We tune the parameters of the simulation 
to achieve the estimated light yield before simulating the 
22 Na source including a flat background contribution to 
account for the events observed above 600 keV from acci- 
dental coincidence of background events in the argon and 
511 keV gamma rays in the Nal detector. We then re-fit 
the data to the simulated spectrum, and the scale factor 
between the integrated charge from both PMTs to the 
energy scale of the simulation provides the energy cali- 
bration; an example of this fit is shown in Fig. [5j While 
the energy resolution is not an explicit parameter of the 
fit, we do smear the photon counts from the simulated 
spectrum with a Gaussian kernel having a standard de- 
viation of X {Npe), as this smearing provides better 
agreement between simulation and data. This additional 
smearing suggests a source of noise in the real experiment 
that is not being accurately modeled in the simulation, 
and a likely source of this noise is the poor gain dispersion 
characteristics of the PMTs. It could also be a property 
of the liquid neon. 

The simulation does not accurately model the data be- 
low about 200 keV, predicting about 10% fewer counts 
than observed down to the low energies where thresh- 
olding effects begin (about 40 keV for the run shown in 
Fig. [SJ. This data excess is seen for a variety of dif- 
ferent gain and trigger threshold settings, and it is cur- 
rently unexplained, although asymmetry effects between 
the PMTs smearing out very low energy events is one 
possibility, as is an additional source of background in 
the detector. Given the discrepancy between data and 
simulation and the fact that that the trigger and gain 
settings were varied throughout the course of the run, 
the energy calibration fit is performed between 200 and 
550 keV. 

There is some tension between the best fit at the 
Compton edge and the 511 keV hump, and the average 
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X^/NDF is 2.5 for all runs with the nominal fit range. 
If the fit is restricted to the range between 350 keV and 
550 keV, the agreement is much better, with an average 
X^/NDF of 1.4. We estimate the systematic uncertainty 
caused by this tension to be 3% by fitting each feature 
independently and taking the difference in the resulting 
energy scale relative to the nominal fit as the uncertainty. 
Performing the fit without the flat background contri- 
bution in the simulation changes the calibration by less 
than 1%. An additional source of uncertainty in this 
calibration are the statistics of fitting the simulated and 
observed spectra, but this source of uncertainty is esti- 
mated to be 1%. 

We cross check the energy calibration obtained from 
the ^^Na source with both higher and lower energy refer- 
ences. For our lower energy reference, we use the ^^Kr™ 
source described in the next section which is the primary 
reason we believe the energy calibration is accurate to 
within our uncertainties. For the high energy reference, 
we look at the Compton edges of the 1.4 MeV and 2.6 
MeV gamma rays produced in the decay of ^"K and ^'^^Tl, 
as there are a number of these events in background data. 
Figure IH] shows background data with the energy scale 
derived from the ^^Na calibration compared to a simula- 
tion of 1.4 MeV and 2.6 MeV gamma rays using the RAT 
simulation. As the purpose of this simulation is solely to 
cross-check the accuracy of the energy scale derived from 
the ^^Na at higher energies, we do not attempt to simu- 
late any other sources of background (i.e. lower energy 
gamma rays in the detector materials, cosmic rays, or 
other backgrounds); therefore, we do not expect the sim- 
ulation to agree with the data below 1 MeV. We find 
evidence for a slight non-linearity in either the energy re- 
sponse of the liquid neon or the PMTs at these energies, 
as the background data do not agree with the simula- 
tion if we assume a linear response of the detector up 
to 2.6 MeV. The fit is improved if the observed energy 
(as calibrated by the ^^Na source) is actually represent- 
ing the following function of the deposited energy in the 
simulation: 



22000 



{Eob 



E, 



dep ^ 



(4) 



It is likely that this dependence is caused by a non- 
linearity of the top PMT, as a comparison using just the 
signal in the bottom PMT requires no energy scaling to 
provide agreement. We also find that the data are best 
fit with an energy resolution of cr = 2AyyEdep. Between 
1 MeV and 2.7 MeV, we find a xVNDF for this fit of 
1.27. 



D. The **'^Kr"' source 

Previously, we reported on use of a ^^Kr"^ generator 
as a source of low energy electronic recoils in liquid argon 
and neon [2^ . This Kr generator became available at the 
end of the experiment described here to provide a second, 
independent energy calibration at 41.5 keV, with a clear 
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FIG. 5: (Color online) The energy calibration is obtained 
by fitting the Compton edge and photoabsorption hump pro- 
duced by 511 keV gamma rays to a simulation of the detector. 
The green lines denote the fit range. We currently have no 
satisfactory explanation for the excess of events between 40 
keV and 200 keV in comparing the data to the simulation. 
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FIG. 6: (Color online) We check the energy calibration ob- 
tained from the '^'^Na source by comparing background data 
to simulations containing ''"K and ^''^Tl components. We 
clearly see the two Compton edges of the 1.4 MeV and 2.6 
MeV gammas produced by these sources. We do not attempt 
to simulate any other backgrounds, so the divergence of data 
and MC below 1000 keV is not surprising. There is evidence 
for a non-linearity in the energy scale at these high energies 
that is likely related to a the response of the top PMT, and 
we correct for this using Eq. [D 



peak appearing in the data as shown in Fig. [7l When 
the energy calibration is derived from the ^^Na spectrum 
as described above, the energy of the ^Sj^j-m pgg^j^ agrees 
with 41.5 keV within statistical uncertainties, suggesting 
both that the energy calibration is accurate and that the 
response of the MicroCLEAN detector is linear between 
9.4 keV and 511 keV (recah that the ^^Ki™ atoms emit 
two electrons at 9.4 keV and 32.1 keV to provide the full 
energy peak). The resolution of our detector at 41.5 kcV 
is 19% (cr/E of the Gaussian fit shown in Fig. [7]). 
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E. Signal yield 

The energy calibration described above uses the total 
signal from both PMTs. Given the absence of a well- 
defined mean single photoelectron in the top PMT, the 
absolute signal yield is estimated from the bottom PMT 
alone by matching the energy calibration to the inte- 
grated charge of the bottom PMT divided by the size 
of the single photoelectron of the bottom PMT. After 
ending the liquid neon run described in this report, we 
filled the detector with liquid argon without warming up 
to room temperature or opening the system; at the higher 
temperatures of liquid argon, we were able to resolve the 
single photoelectron response of both PMTs. From these 
measurements, we determined that the top tube observes 
(6 ± 3)% more signal than the bottom tube [ll], due 
to some combination of photon collection efficiency and 
photocathode quantum efficiency. Therefore, we extrap- 
olate the total light yield in the liquid neon to include 
both tubes based on the observed PMT efficiencies. Us- 
ing this method, we determine a photoelectron yield of 
(3.0 ± 0.3) photoelectrons/keVee from the ^^Kr"^ data. 
The total uncertainty includes an estimate of the uncer- 
tainty in the relative efficiency of the two tubes intro- 
duced by changing the temperature from 85 K to 25 K, 
but it is dominated by the single photoelectron response 
of the bottom PMT. 

We next use the ^^Kr™ data to validate the ^^Na cali- 
bration method. Neglecting the systematic uncertainties 
from the PMT photoelectron response and efficiency, as 
they apply equally to both methods for determining the 
light yield, the standard ^^Na energy calibration taken 
at the same time as the ^^Ki"^ data returns a light yield 
of (3.00 ±0.02) photoelectrons/keVee, in agreement with 
the 83j^-|.m calibration. This agreement provides confi- 
dence that we can trust the ^^Na calibration for the data 
taken before the the ^^Ki"^ source became available. 

As will be discussed in Sec. IIIFi the signal yield was 
not stable throughout the course of the run; the largest 
signal yield observed in this study was (3.5 ± 0.4) pho- 
toelectrons/keVee, a factor of about 3.5 larger than we 
achieved in our previous measurement (0.9 photoelec- 
trons/keVee We attribute the increase to immersing 
the PMTs in the liquid instead of viewing through pres- 
sure windows, improved photocathode coverage, and a 
change in the thickness of the wavelength shifter. 
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FIG. 7: (Color online) Energy spectrum of a background sub- 
tracted *^Kr™ run in neon. As discussed in the text, the 
observed light yield for these data is (3.0 ± 0.3) photoelec- 
trons/keVee. The resolution is 19% (c/E of the Gaussian fit 
shown in red) at 41.5 keV. 



gaged the circulation pump and flowed neon through the 
charcoal trap shown in the schematic of Fig. [2] to remove 
impurities. The light yield subsequently increased, reach- 
ing a peak at (3.5 ± 0.4) photoelectrons/keVee as shown 
in Fig. [HI providing evidence that impurities were being 
contained by the charcoal. Because we could not purge 
the internal charcoal trap while running the experiment, 
it is not clear whether all the impurities in the neon were 
removed or whether the trap became saturated. These 
data were taken for a liquid temperature of 28.7 K. 
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F. Detector stability 

The presence of impurities can reduce the amount of 
light detected both by quenching excimer states non- 
radiatively or by absorbing the VUV scintillation light, 
and both processes have been studied in liquid argon |29l - 
[3l| . We monitored both the light yield and the triplet 
lifetime over the course of the experiment to gauge the 
effect of impurities. Initially, the signal yield of the de- 
tector was only 1.9 photoelectrons/keVee. We then en- 



FIG. 8: The light yield in the detector as a function of time 
with the circulation pump to the charcoal trap engaged for 
a liquid temperature of 28.7 K as determined by the ^^Na 
source. The light yield reached a maximum at (3.5 ± 0.4) 
photoelectrons/keVee. The error bars are statistical only and 
do not include the uncertainty introduced by the single pho- 
toelectron response of the bottom PMT. 

Both the light yield and the triplet lifetime varied with 
time. As will be discussed in Sec. |Vl both of these pa- 
rameters are strongly and unexpectedly affected by the 
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temperature and pressure of the liquid neon. To sta- 
bilize the temperature and pressure, we disengaged the 
circulation pump for the remainder of the run, poten- 
tially allowing contaminants to build up in the neon and 
decreasing the observed light yield in the 83j^^.m jg^^g^ 
taken at the very end of the experiment to the (3.0 ±0.3) 
photoelectrons/keVee reported above. 

To account for the temperature effects, the PSD and 
£eff data were collected at three different temperatures, 
26.7 K, 27.8 K and 28.7 K, where each data set was 
taken over a roughly two week period. We measured 
the light yield approximately every other day during 
each run to determine the stability of the detector for 
each temperature set point. The observed signal yields 
were 2.74 ± 0.03, 3.15 ± 0.08 and 3.13 ± 0.03 photoelec- 
trons/keVee at 26.7 K, 27.8 K and 28.7 K respectively, 
where the error bars represent the root-mean-square er- 
ror of all the calibration runs taken at the operating tem- 
perature. We therefore conclude that the detector signal 
yield was stable to better than 3% for each temperature 
setting. 
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FIG. 9; (Color online) Mean prompt fraction as a function 
of energy for electronic and nuclear recoils at three dilferent 
temperatures. The error bars include both statistical and 
systematic uncertainties and are generally dominated by a 
2.5% uncertainty stemming from dilferences in the two PMTs. 



III. PULSE SHAPE DISCRIMINATION 



We define the prompt fraction /„ as 



(5) 



where V{t) is the vohage trace from the PMT, ^ = 220 
ns is an integration period determined to optimize the 
PSD Ql,T,^tQ- 50 ns, Tf = to + 44 /is, and to is the 
trigger time. As mentioned in the previous section, we 
collected data at 26.7 K, 27.8 K and 28.8 K. In energy 
bins of 5 keVee, we fit a Gaussian function to the empiri- 
cal distributions to determine the mean prompt fraction, 
/p, for nuclear and electronic recoils, and Figure [3] shows 
fp as a function of energy for both classes of events. The 
error bars include both statistical uncertainties and a sys- 
tematic uncertainty of 2.5% stemming from differences in 
the two PMTs, with the systematic error generally dom- 
inating; the total uncertainty is smaller than the size of 
the markers for the electronic recoil data. There is also 
a small systematic error associated with fitting an asym- 
metric distribution with a symmetric Gaussian function 
that we do not include in this analysis. The nuclear re- 
coil data were acquired in coincidence with a BC501A 
organic scintillator module as described in Sec. Ill Bl Two 
effects may be observed: first, the mean prompt fractions 
for nuclear and electronic recoils converge at low ener- 
gies. This effect has been observed in liquid argon [9l.[l0l| 
and is likely due to the fact that dE/dx for nuclear and 
electronic recoils also converge as the energy decreases. 
Also, the mean prompt fraction decreases with increasing 
temperature and pressure for both electronic and nuclear 
recoils at all energies. This will be discussed further in 
Sec.lVl 



We define the electronic recoil contamination (ERG) 
as the probability of mistaking an electronic recoil for 
a nuclear recoil given a particular level of nuclear re- 
coil acceptance. We estimate the ERG as the number of 
tagged electronic recoil events with fp > /p, nuclear, set- 
ting the nuclear recoil acceptance level to approximately 
50%. The observed ERG at each temperature is shown 
in Fig. [TOl Our observed ERG is about a factor of 5- 
10 better than that observed in Q. Theoretically, we 
might have expected a stronger improvement given that 
we collect at least three times as much light. The smaller 
than expected improvement is in part attributable to the 
use of a /p, nuclear that chauges with energy, instead of 
the flat cut value assumed previously in [6|. As shown 
in Fig. [9l the mean prompt fractions for electronic and 
nuclear recoils converge at lower energies, reducing the 
discrimination power. 
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FIG. 10: (Color online) ERG versus energy for the three data 
sets. The x-axis is given in relevant units for dark matter 
searches, keVr, by assuming jCeH = 0.24 (see Sec. IIV|). 
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A second cause of weak PSD is the electronic noise in 
the experiment. We fit the data using a statistical model 
for fp described in 0, , and Figure [TT] shows the PSD 
for the 26.7 K data with a fit to this model. In Q, we 
described a method for estimating the size of the ad- 
ditional noise for comparison between detectors, model- 
ing the prompt and late noise as a constant multiplied 
by the mean prompt and late signals: dp = Cp^p and 
af — Cifii. As discussed in Q, we are wary of interpret- 
ing these fit parameters as a literal measurement of the 
noise; however, we do believe they can be used as points 
of comparison. With that method, the noise parameters 
in the liquid neon are a factor of two to three larger than 
those observed in argon. The probable cause for such a 
large amount of noise is the poor gain characteristics of 
the PMTs and the use of an additional amplifier phys- 
ically far away from the PMT itself, in addition to the 
factor of four longer integration period required for liq- 
uid neon scintillation. The solid line in Fig. [TT] shows the 
statistical model prediction in the ideal case where the 
additional noise is set to zero. We note that the effect of 
the additional noise is most relevant for energies above 
100 keVr. 




Energy (keVr) 



FIG. 11: PSD for the 26.7 K data set individually, along with 
a fit to the statistical model. The solid line shows the model 
prediction with additional noise terms set to zero. We note 
that the effect of the additional noise is most relevant for 
energies above 100 keVr. 

To test the hypothesis that the poor PMT response 
negatively affected the PSD results, the signals in each 
PMT were analyzed separately. As expected given the 
single photoelectron distributions, the bottom PMT pro- 
vided as much as an order of magnitude better PSD com- 
pared to the top PMT above 150 keVr, dropping to a fac- 
tor of about 2 at lower energies, although the combina- 
tion of the two PMTs still provides an improvement over 
just the bottom PMT alone. To achieve better overall 
PSD, more effort will be needed to improve the response 
of the PMTs at liquid neon temperatures; as the two 
PMTs exhibited different behaviors when cold, all PMTs 
intended for use in a liquid neon detector may need to 



be cryogenically tested. Another potential area of im- 
provement would be a more sophisticated discrimination 
statistic. Given the complicated nature of scintillation in 
liquid neon (see Sec. |Vl), a PSD analysis using more of 
the timing information along the lines of the multi-bin 
method described in Q might prove rewarding. 



IV. NUCLEAR RECOIL SCINTILLATION 
EFFICIENCY 

As mentioned in the introduction, it is well known that 
nuclear recoils produce less scintillation light in liquid 
noble gases than electronic recoils do. There are several 
processes related to the linear energy transfer (LET) of 
the incoming radiation that are known to contribute to 
the lower scintillation efficiency of nuclear recoils. The 
first of these is the Lindhard effect (s^l , where some of the 
energy deposited goes into heat instead of the creation 
of molecules. The Lindhard effect has been shown to be 
sufficient to explain quenching in germanium [s^, but 
it is inadequate to explain CeS for the noble gases. A 
second mechanism is known as bi-excitonic or Hitachi 
quenching [34|. whereby 

Ne* + Ne* ^ Nc + Ne+ + e" . (6) 

When such a reaction occurs, the potential for creating 
two photons from the two original exciton states has been 
reduced to one photon via recombination of the single re- 
sulting ion. This reaction is more likely to occur in high 
density tracks like those produced by nuclear recoils. Re- 
cently, a model has been proposed to phenomenologically 
tie this quenching to the stopping power [l^, [sl] . Lastly, 
some fraction of ionized electrons never recombine; if 
these "escape electrons" are more likely for nuclear re- 
coils than electronic recoils, this process would contribute 
to a smaller CcS- This effect has been observed to lower 
£cS in liquid xenon [l^ |3^, [s^l . 

The experimental setup for measuring the nuclear re- 
coil scintillation efficiency has been described in [l^ in 
the context of a measurement in liquid argon. One 
salient difference between the current analysis and the 
argon analysis is the absence of MC simulations of the 
neutron scattering. There are few experimental data 
on neon-neutron scattering, with Geant4 switching dis- 
cretely from data to a model for neutrons of energy < 20 
MeV. There are efforts underway to improve this situa- 
tion for future studies [sst . 

Instead of fitting MC spectra to the data to determine 
the scintillation efficiency, a simple Gaussian function is 
fit to the data to determine the mean of the observed 
peak in units of keVee. Given the non- Gaussian nature 
of the single photoelectron distribution, we expect some 
skewing of the observed spectra towards high energies. In 
a simple MC of a similar situation in a larger detector, 
the contribution of this skew is « 5% for the lowest en- 
ergy point, below the level of our dominant uncertainties. 
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Hence, our scientific judgment is that the uncertainty in- 
troduced by this skew does not contribute significantly to 
our final uncertainty estimate. The final expression for 
CcS is the observed mean of the signal yield divided by 
the expected energy of the recoil as predicted by Eq. [3l 

Data were taken at nine different scattering angles cor- 
responding to energies ranging from 28.9 keVr to 368.7 
keVr. Each point, with the exception of the point at 
178 keVr, was taken at the three different operating tem- 
peratures. The results at each temperature are consis- 
tent with each other to within the uncertainties of the 
measurement, so for the final determination of CeS, all 
the data from each temperature are combined into one 
data set. Figure [T^ shows the final Gaussian fits at each 
energy, and Figure [T3] shows the observed nuclear recoil 
scintillation efficiency, along with the result from the pre- 
vious measurement at 387 keVr Q. To obtain Ccs we 
assume that the signal yield from electronic recoils scales 
linearly with energy (consistent with our observations be- 
tween 40 keVee and 511 keVee) and compare to the cor- 
responding signal yield at each nuclear recoil energy. The 
mean value of Ccs thus obtained is 0.24. 

Because we expect there to be some contamination 
from multiple-scattering and gamma backgrounds into 
the distribution, as observed in argon data [1^ , we limit 
the fit range in the region of the peak in the distribution; 
the fit range is iteratively chosen to encompass plus and 
minus one sigma from the central value representing our 
estimate of the mean energy of the peak. The uncertainty 
associated with this choice of fit range will be discussed 
in Sec. IIVBI The average x^/NDF for the fits shown in 
Fig.iniis 1.07; the largest value of xVNDF is 2.5 for the 
132.2 keVr peak, with all others below 1.3. 



A. Data selection cuts 

One data selection cut requires that the two PMTs 
trigger within 40 ns of each other and that the asymmetry 
parameter is within one sigma of the mean asymmetry 
value for the run. To select neutron-like events, a PSD 
cut was applied to the organic scintillator data as shown 
in Fig. [M] We also apply a TOF cut between an event in 
the neon and the tagging event in the organic scintillator. 
A final cut is applied based on the PSD described in 
the previous section; this is a very loose cut, eliminating 
events farther than two sigma from the measured mean 
value of fp, for nuclear recoils at 26.7 K. 

The TOF cut appears to be quite effective in selecting 
single-scatter nuclear recoil events in liquid neon. For 
all runs, neutron events formed a clear peak in the TOF 
spectrum. As single-scattering neutrons are more likely 
to be found at the very beginning of the TOF peak, 
the TOF cut was located at the front side of the peak. 
First, the total number of counts in the neutron peak 
was found, and then the TOF distribution was scanned 
until the number of counts in a 4-ns-wide bin exceeded 
3% of the total integral of the neutron peak; this point 



set the left boundary of the TOF window. The width of 
the TOF acceptance window was chosen to be 10 ns. An 
example of the TOF cut is shown in Fig. [Tn]for 368.7 keV 
events. 

Figures [16] and [17] show the change in the observed 
energy spectra produced by sliding the location of the 
TOF cut window to earlier and later times for both the 
lowest and highest energy runs. If the window is moved 
10 ns earlier, very few events pass the cut and no peak 
in energy is visible. At the standard cut location, clear 
peaks in the energy spectra for both the low energy and 
high energy points are visible. As the window moves to 
accept events with longer times-of-fiight, the energy peak 
loses definition due to the inclusion of multiple-scattering 
backgrounds, becoming an almost fiat background when 
the TOF cut is moved to the center of the TOF peak (at 
20 ns beyond the standard cut). These results provide 
confidence that the standard TOF cut at the front side 
of the TOF peak is selecting single-scattering events and 
eliminating multiple-scattering backgrounds. 



B. Error analysis 

There are a number of sources of uncertainty in this 
measurement of CeS relating to both experimental pa- 
rameters as well as data selection cuts. There is a general 
uncertainty of 3% from the energy calibration described 
above. Studies of the analysis cuts were performed in- 
dependently to assess the size of the systematic error as- 
sociated with each. In each case, we took the difference 
between the value of Ces determined by the standard 
analysis and that determined with the changed cut value 
as our estimate of the systematic error associated with 
that cut. The size of the coincidence window between the 
trigger times of the two PMTs was narrowed to 20 ns, and 
the cut was dropped entirely to study its effect. Both the 
individual trigger cuts and the asymmetry cut were ex- 
panded to encompass the region within two sigma of the 
mean for the entire run. The range in pulse area of the 
organic scintillator PSD cut was first narrowed and then 
scanned across the entire distribution of nuclear recoils 
(the lower band in Fig.[T4]). To assess the uncertainty as- 
sociated with the TOF cut, the width of the TOF window 
was both narrowed to 6 ns and expanded to 16 ns, with- 
out changing the location; as already discussed, moving 
the TOF cut window away from the front of the TOF 
peak eliminated any observable peak, so the location of 
the cut was not included as a systematic uncertainty. Fi- 
nally, the analysis was repeated with no PSD cut applied 
to the neon data. To estimate the size of the uncertainty 
associated with the fit range in Fig. [T^ the range was 
expanded to include the region within two sigma of the 
central value. The uncertainty associated with the elec- 
tron equivalent energy calibration was determined to be 
2%. A source of uncorrelated error is uncertainty in the 
location of the organic scintillator affecting the predicted 
value of the nuclear recoil energy. An uncertainty of 0.64 
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FIG. 12: (Color online) Fits for the nuclear recoil scintillation efficiency in liquid neon for all energies. All data sets include 
results from each temperature setting, except for the data at 178 keVr that were only taken at 26.7 K. 
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FIG. 13: (Color online) The observed nuclear recoil scintilla- 
tion efficiency versus nuclear recoil energy in neon, along with 
the Lindhard+Birks model described in the text. 



FIG. 14: (Color online) The PSD cut in the organic scintilla- 
tor to select neutron events. Only events within the pentagon 
are included in the data set. 
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FIG. 15: The TOP cut window was 10 ns, and only events 
at the beginnning of the TOF peak were selected, as these 
events are more likely to be single-scatter events. 

cm in the location of the organic scintillator around the 
detector corresponds to an angular uncertainty of 1.3° 
at each position. A final source of uncertainty is the 
trigger efficiency. Data for the three lowest energy posi- 
tions were taken at different hardware trigger thresholds 
such that the difference between the lowest and highest 
threshold was larger than the average height of a single 
photoelectron in each PMT. 

The statistical uncertainties in the fit results of Fig. [12] 
were generally below the level of the other uncertainties 
in the measurement. The main source of uncertainty 
for the lower energies found in the experiment was the 
choice of fit range, producing an uncertainty of as large 
as 15% in the 28.9 keVr data set. In all cases, the effect 
of increasing the fit range was to increase the observed 
value of iZeff. The results from changing the hardware 
threshold were consistent with each other at each of the 
lowest three energy points. At the highest energies, the 
dominant source of uncertainty was the width of the TOF 
window. All sources of uncertainty were assumed to be 
independent and combined in quadrature to produce the 
final uncertainty values shown in Fig. [131 Table [ij lists 
the central values and the uncertainties associated with 
each identified source of error. 



C. Discussion and the Lindhard+Birks Mode 

The first point to note about our results for Ccs in 
liquid neon is an upturn at low energies similar to one 
observed in our previous measurement of CcS in liquid 
argon [l5j, although the uncertainties on the lowest en- 
ergy points are also the most significant. Both hardware 
and software trigger effects were examined in attempts to 
explain this upturn, with neither producing the observed 
result. As with the argon, we do not know whether this 
is a physically real effect. 

In the introduction, two main processes were listed as 



causes of the lower scintillation efficiency of nuclear re- 
coils in liquid noble gases. The first was the Lindhard ef- 
fect, where some energy goes into heat instead of the cre- 
ation of molecules. A second effect was Hitachi quench- 
ing, where two excited atomic states interact to produce 
one ground state atom and one ion, reducing the total 
number of Ne2 molecules available to produce scintilla- 
tion photons. Mei et al. have developed a phe nomeno- 
logical model to account for these two effects [331 > s-'^d 
apply that model here to the current results drawing val- 
ues for the stopping power from the tables of SRIM [3^ 
("SRIM" in Fig. [T3|). We also use stopping powers ob- 
tained from Mei et al. from a second model for stopping 
power as described in ^ ("DM" in Fig.[l3|). 

The result of fitting the Lindhard+Birks model to our 
data is shown in Fig. [TSj The Birks parameters, fcS, for 
the SRIM and DM stopping powers are determined to 
be 1.99 X 10-3 and 1.73 x lO'^ MeV'^ g cm'^, respec- 
tively. This can be compared to the value of 1.12 x 10"'^ 
MeV""'^ g cm~^ determined in [s^ from the single data 
point marked as "PicoCLEAN result" in Fig. [T3 The 
kink that appears around 40 keVr in the SRIM curve is 
potentially mirrored in the upturn at low energies but is 
not present in the DM model. In developing the Lind- 
hard+Birks model for argon, Mei et al. found a similar 
kink from the SRIM data that was not present in alter- 
native models [35^ . 



V. SCINTILLATION TIME DEPENDENCE 

The established model for scintillation in liquefied no- 
ble gases consists of two exponentially decaying com- 
ponents representing the contribution from singlet and 
triplet molecules. In our measurements of scintillation 
in liquid neon, although fast and slow exponential be- 
havior is observed, neither a two-component nor a three- 
component model is sufficient to describe the time de- 
pendence. In addition, the observed time constants and 
intensities of each component vary significantly with the 
temperature and pressure of the liquid. In an attempt to 
understand the various components of liquid neon scin- 
tillation and the temperature dependence, we examine 
average traces from electronic recoil events with energies 
between 300 keVee and 400 keVee for nine different pres- 
sure and temperature settings. All data are derived from 
tagged ^^Na events; the trigger times of each event in the 
data sets are aligned and the events are summed together 
to produce the average traces. 

Figure [TH| shows an example trace for data taken at 
27.4 K, along with fits to a two- and three-exponential 
model. While the three-exponential model improves on 
the quality of the overall fit, it still does not accurately 
model the data between about 50 ns and 1 ^s. We there- 
fore adopt a four-exponential mixture model |40j. The 
exponential decay constants are denoted Tj, and the rel- 
ative weights, which are non- negative and sum to 1, are 
denoted li. The mixture model for the probability den- 
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FIG. 16: The observed energy spectrum for 368.7 keVr events as the location of the TOF cut (as defined by TOFl and TOFh) is 

moved by ±10 ns and +20 ns. The nominal cut around the front side of the TOF peak is the case where TOFl < TOF < TOFh. 
The degradation in the observed peak as the TOF cut moves away from standard location is clear. 
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FIG. 17: The observed energy spectrum for 28.9 keVr events as the location of the TOP cut (as defined by TOFl and TOFh) is 
moved by ±10 ns and +20 ns. The nominal cut around the front side of the TOF peak is the case where TOFl < TOF < TOFh. 
The degradation in the observed peak as the TOF cut moves away from standard location is clear. 
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TABLE L The observed L^b parameters as a function of energy along with the estimated uncertainties for each of the sources of 
uncertainty listed in the text. The fit range column refers to expanding the fit range out to two sigma. The hardware threshold 
was varied for only the three lowest energy points. The "OSc cut" column refers to the location of the data cut used on the 
organic scintillator neutron detector. The "TOF size" refers to the size of the time-of- flight cut, which was both expanded to 16 
ns and shortened to 6 ns, with both contributions included in that column. The "Temp." column lists the standard deviation 
of the three temperature set points. All uncertainties that did not contribute at the hundredths level have been listed as zero. 
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FIG. 18: (Color online) The average traces for 27.4 K data, 
along with fits to two-component and three-component expo- 
nential models. Neither model accurately describes the inter- 
mediate time regime. 



sity function of the creation times of scintillation photons 
created by an energy deposition at t = is 



/(t) = E^e-*/-. 



(7) 



We model the expected voltage trace produced by the 
PMTs as a convolution of an experimentally determined 
impulse response function of the PMT, h{t), and the PDF 
model: 



{Vit)) (X / h{t-s)f{s~to)ds, 



(8) 



where Iq is the time of the energy deposit. We perform 



the fit by minimizing the Matusita distance|4l| between 
the model and the data. 

Example fits to the 26.7 K and 28.8 K data are shown 
in Fig. [T21 We associate ti and T2 with the longest 
and shortest time components of the scintillation, respec- 
tively, and they are shown in Fig. [20] as a function of 
neon temperature. The short time constant, T2, should 
not be strictly associated with the decay of the singlet 
state, as the fast component is also affected by the tim- 
ing characteristics of the TPB wavelength shifter and the 
PMTs. The remaining scintillation time constants, rep- 
resenting intermediate time scales, are shown in Fig. 1211 
The weights of each component as a function of temper- 
ature are shown in Fig. [22l and all fit parameters are 
reported in Table [TTl 

We estimate one sigma random uncertainties using a 
bootstrap resampling scheme [4^ . We determine approx- 
imate estimates of the systematic uncertainties associ- 
ated with choice of fit window by varying the end time of 
the integral in Eq. [5]between 44 and 52 fj,s. We estimate 
the systematic uncertainty associated with our measure- 
ment of the single photoelectron response of the PMTs by 
using several different estimates of h{t) taken from differ- 
ent data sets throughout the run. We also perform the fit 
for each PMT individually to estimate the uncertainties 
associated with the different PMTs. In Figs. BOll^ the 
error bars represent the combination in quadrature of the 
statistical and systematic uncertainties. In many cases, 
the error bars are smaller than the size of the marker. As 
a further consistency check, the analysis was repeated for 
events with energy between 50 and 150 keVee, and all the 
trends discussed below were observed. 

We make a few general observations regarding the data 
before continuing to a discussion of the physics. First, the 
observed time constants decrease with increasing tem- 
perature and pressure (Figs. [^OlET]) . Second, the rela- 
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FIG. 19: (Color online) Examples of the four-component 
model fit to 26.7 K and 28.8 K data. 
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FIG. 21; (Color online) The fitted intermediate time con- 
stants, T3 and T4, for neon scintillation as a function of neon 
temperature. The error bars represent the combined esti- 
mated statistical and systematic uncertainties, derived as de- 
scribed in the text. 
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FIG. 22: (Color online) The fitted weights of the four- 
component model for neon scintillation as a function of neon 
temperature. The error bars represent the combined esti- 
mated statistical and systematic uncertainties, derived as de- 
scribed in the text. 



FIG. 20: (Color online) The fitted longest and shortest time 
constants, n and T2, for neon scintillation as a function of 
neon temperature. The error bars represent the combined 
estimated statistical and systematic uncertainties, derived as 
described in the text. 



live intensity of the long-lived component increases with 
temperature, primarily at the expense of the interme- 
diate components with a small contribution taken from 
the short-lived part (Fig. [22]). The total signal yield also 
changes with temperature. Figure [231 shows the contri- 
bution to the total signal yield by each of the four ex- 
ponential decay processes in the model (A,; = x Ajoj 
where Xtot is the total signal yield) at four temperatures 
taken within three days to minimize the effect of chang- 



ing detector conditions; the total signal yield increases 
with increasing temperature, suggesting that not all of 
the increase in Ii is due to the decreases in the other 
components. 

There are two final observations that can be made from 
the data. First, the long time constant is affected by the 
presence of impurities, as the long time constant at con- 
stant temperature increased during the purification cycle 
described in Sec. Ill Fl Second, we perform a similar anal- 
ysis on untagged nuclear recoil data at two different tem- 
peratures (unfortunately, the tagged nuclear recoil data 
are limited by statistics) and present the results at the 
bottom of Table |TT1 The parameter trends between the 
two nuclear recoil data sets are consistent with the elec- 
tronic recoil data, although the long time constant for 
nuclear recoil events appears to be 1 — 2 /xs shorter than 
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T p Ti (Tti T2 cr t2 ra (Tra r4 (7^4 /i tj/j J2 cr/^ J3 cr/., J4 crj^ 
(K) (mbar) (/is) (ns) (/is) (/is) 

Electronic recoils 

26.7 816 18.37 1.13 17.93 0.70 3.00 0.12 0.31 0.03 0.56 0.02 0.16 0.01 0.23 0.02 0.049 0.002 

27.4 1048 18.34 1.12 17.05 0.75 2.72 0.18 0.26 0.03 0.66 0.02 0.13 0.01 0.17 0.02 0.039 0.003 
27.7 1160 16.14 1.10 16.22 1.02 2.23 0.43 0.20 0.04 0.73 0.03 0.11 0.01 0.12 0.02 0.032 0.004 
28.0 1281 16.40 1.05 16.50 0.86 2.32 0.36 0.22 0.04 0.71 0.03 0.12 0.01 0.13 0.02 0.034 0.004 
28.3 1380 13.77 0.15 15.72 0.74 1.46 0.37 0.15 0.02 0.81 0.02 0.097 0.005 0.07 0.01 0.025 0.004 

28.5 1467 13.60 0.25 14.97 0.78 1.48 0.29 0.14 0.02 0.81 0.02 0.10 0.01 0.07 0.01 0.025 0.003 

28.6 1520 12.33 0.26 14.88 0.51 0.89 0.12 0.10 0.01 0.85 0.01 0.089 0.004 0.04 0.01 0.020 0.002 

28.7 1575 12.09 0.38 13.77 0.67 0.75 0.08 0.09 0.01 0.86 0.01 0.08 0.01 0.037 0.004 0.020 0.002 
28.9 1662 12.07 0.39 13.95 0.73 0.73 0.08 0.09 0.01 0.86 0.01 0.08 0.01 0.035 0.004 0.019 0.002 



Nuclear recoils 


27.7 


1160 


13.80 


0.75 


18.38 


0.88 2.29 0.34 0.31 0.05 0.43 


0.01 


0.35 


0.01 


0.142 


0.005 


0.08 


0.01 


28.3 


1380 


10.99 


0.51 


17.49 


1.02 1.26 0.31 0.175 0.04 0.52 


0.02 


0.32 


0.01 


0.105 


0.004 


0.05 


0.01 



TABLE II: Estimated time dependence parameters of the four-exponential model for scintillation in liquid neon. The long and 
short components are ti and T2 respectively, and the intermediate components are ra and T4,. h represents the relative weight 
of each component. 
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FIG. 23: (Color online) The total signal yield plotted with 
the intensities of the long, short, intermediate and residual 
components as a function of temperature multiplied by the 
total signal yield of four runs taken over the course of three 
days. 

the corresponding time constants for electronic recoils. 

A. Time dependence discussion 

Besides the very prompt light produced by the de- 
cay of singlet states, there are three additional, distinct 
scintillation components. The triplet lifetime in liquid 
neon has previously been measured to be 15.4, 2.9 and 
3.9 /Its [Ijli^lll], with the large discrepancy between the 
first and second two measurements attributed to impu- 
rities in the liquid in the latter two. In solid neon, the 
long exponential time constant has been measured to be 
3.3 and 3.9 /j,s [i^ HJ]. A very slow tail in the solid has 



also been observed with a lifetime of 560 us, attributed 
to excitations of the 3sf2P2 atomic state [8|; in one mea- 
surement of the liquid, an excited atomic triplet state was 
observed to decay non-exponentially, disappearing within 
several microseconds fiSi] . Measurements in gaseous neon 
suggest that long-lived ^Pi and '^Pi atomic states can re- 
act via three-body collisions to form radiative molecular 
states [4^ ; the lifetime of the long-lived molecular state 
in gas, presumably equivalent to the triplet state, was 
measured to be 6.6 /iS [i^. Lastly, as the temperature of 
the gas is reduced, the emission due to molecular states 
disappears, leaving only the atomic spectra behind 8]. 

Before further discussing the current neon data, we 
compare the results to the situation in argon and he- 
lium. First, the data presented in [9|] for scintillation in 
liquid argon were taken at many different temperatures 
covering a larger absolute temperature range than that 
examined in liquid neon, and there was no change in the 
observed triplet lifetime or signal yield in argon. 

In liquid helium, the triplet lifetime is 13 s [l^l- Two 
other characteristic time scales were observed in scintil- 
lation in liquid helium. First, metastable He(2'^S) atoms 
created by ionizing radiation have been observed to last 
for 8000 s in vacuum |4§|; in liquid, these states only 
last for 15 /xs, and their disappearance is accompanied by 
the appearance of vibrationally excited triplet molecules, 
implying the creation of molecules by collisions of the 
excited atomic states [50| . Second, a component has 
been observed over a finite time window and attributed 
to collisions of long-lived triplet states leading to the pro- 
duction of the shorter-lived singlet states which promptly 
decay [5l| . King and Voltz developed a model to describe 
slow emission light produced by scintillators via triplet- 
triplet collisions, neglecting the effects of the radiative de- 
cay of triplet molecules [52]. The KV model was adapted 
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by McKinsey et al. to describe the component in he- 
lium over a finite time window, although the exact form 
of t he p ower law depends on the assumed track geome- 
try 

Finally, an important similarity between liquid helium 
and liquid neon that sets them apart from the heavier no- 
ble gases is that free electrons form bubbles or "localized" 
states in both hquids [53l - [57| . These bubble states occur 
because the electron-atom repulsion due to Pauli exclu- 
sion is strong relative to the weak polarization of helium 
and neon atoms, as evidenced by a positive scatte ring 
length for low energy electrons in helium and neon |58l - 
l60j . In the heavier noble gases, the polarization term be- 
comes stronger and overwhelms the repulsive pseudopo- 
tential. While electron bubbles would influence recombi- 
nation rates, affecting the timing of scintillation on short 
time scales as well as the number density of singlet and 
triplet molecules, it is unlikely that they are important 
for the late time behavior. However , in liqu id helium, Hcj 
molecules also form bubble states |6ll - l63j . Similar stud- 
ies and calculations have not been performed for liquid 
neon, but it is certainly possible that Ne2 could also form 
bubble states, affecting the long-lived triplet molecules 
and providing another difference between scintillation in 
liquid neon and liquid argon. 

The comparison to liquid helium seems to provide 
the most illumination, particularly as free electrons form 
bubbles in both liquids. As in the current measurements, 
four time components have been observed in liquid he- 
lium scintillation: a prompt decay from singlet molecules, 
a slow exponential decay from triplet molecules, an inter- 
mediate exponential decay from atomic states that can 
also react to form molecules, and a component from 
triplet molecules reacting to produce singlet molecules. 

In neon, we likewise attribute the slowest and fastest 
(ti and T2) components to singlet and triplet molecules. 
The number of promptly created singlet molecules de- 
creases with increasing temperature and pressure, while 
the number of triplet molecules greatly increases. We 
know that there is a connection between the recoil 
track density and the production of singlet and triplet 
molecules, and a possible contributor to the changing in- 
tensities of each component is the changing density of ex- 
citations with temperature. From measurements of the 
electron mobility and the Einstein relation, the diffusion 
constant for localized electrons increases by 30% between 
26.7 and 28.8 K [s^. Recalling that one channel for the 
production of metastable molecules is ion-electron recom- 
bination, in a simplified picture one expects a larger diffu- 
sion constant for electrons to result in a greater maximal 
separation between ions and recombining electrons. This 
could produce a greater proportion of triplet molecules 
given the necessity of an electron spin flip, a mechanism 
that has been proposed to explain the behavior of xenon 
scintillation under an electric field [l3|. However, the in- 
crease in Ii with increasing temperature is much greater 
than the corresponding decrease in I2, requiring some 
additional explanation that may be found in the inter- 



mediate components. 

Continuing to examine the triplet states, we must 
also conclude that the apparent lifetime of the triplet 
molecules changes as a function of the temperature or 
pressure of the liquid environment. Impurities in liquid 
noble gases can shorten the apparent triplet lifetime via 
absorption or non-radiative collisions. However, in the 
data of Table [TTl the shorter lifetime observed at higher 
temperature and pressure is not accompanied by a de- 
crease in Ii; in fact, the contribution of the triplet com- 
ponent greatly increases with decreasing ti, contrary to 
what one would expect if decreasing ti was caused by 
impurities. 

One possible mechanism that could explain the chang- 
ing triplet lifetime without an accompanying reduction in 
Ii is the presence of molecular bubble states analogous 
to those found in liquid helium; one could easily imag- 
ine that the radiative lifetime of the molecule changes 
depending on whether it resides within a bubble or not. 
If Ne2 molecules are weakly bound in bubbles at these 
temperatures, then some fraction of molecules could be 
localized in bubbles with an additional fraction freely ex- 
isting in the liquid, leading to the following model for the 
decay rate of triplet molecules, d{Nt)/dt: 



dt 



n 



n 



Here, f{p,T) is a function of temperature and pressure 
representing the fraction of molecules contained in bub- 
bles, Tb is the lifetime of molecules in bubbles and ti is the 
lifetime of molecules outside of bubbles. Although such 
a model could account for the observed behavior, it is 
purely speculative given our lack of theoretical guidance 
for whether neon molecules truly do form stable bubbles. 

There are two intermediate components, one with a 
lifetime of order 100 ns and the other with a lifetime of 
order 1 /xs. We can again draw a parallel to the situa- 
tion in liquid helium and attribute some of this compo- 
nent to collisions between triplet molecules seeding the 
creation of singlet molecules that immediately decay as 
in the King and Voltz theory. However, the KV theory 
cannot entirely account for the intermediate exponential 
behavior, even if one modifies the KV model to include 
the decay of the triplet state. Therefore, we attribute 
the longer-lived intermediate exponential component to 
either excited atomic states producing singlet molecules 
that immediately decay or to the decay of the long-lived 
atomic states themselves. In the first case, the observed 
lifetime T3 is equivalent to the reaction rate of excited 
atomic states with ground-state Ne atoms; if this collision 
rate increases with temperature or pressure (perhaps be- 
cause of increasing diffusion or bubble fluctuations), the 
observed lifetime would decrease. In the latter case, T3 
would be some combination of the true radiative lifetime 
of the atomic excitation and a quenching reaction. In 
contrast to the very long-lived exponential channel, the 
decreasing ra is accompanied by a decreasing , suggest- 
ing the existence of a non-radiative channel. Taking the 
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argument one step further, if we combine the two sce- 
narios with the modification that three-body colHsions 
of excited atomic states create triplet molecules, the non- 
radiative channel could result in the production of triplet 
molecules and the observed increase in Ii . This interpre- 
tation could be consistent with the disappearance of a 
triplet atomic spectrum over several microseconds as ob- 
served by Suemoto and Kanzaki 43] . 

In the future, understanding the effects of an applied 
electric field on liquid neon scintillation might illuminate 
some of the temperature/pressure dependencies observed 
in this work. A second study that would be useful in 
assessing whether molecules reside in bubble states 
would be further laser spectroscopy; in liquid helium, 
molecular bubbles were first discovered in the absence of 
wavelength shifts for rovibrational structure in emission 
spectra of He2 between gaseous and liquid helium. These 
studies could be accompanied by theoretical work on the 
potential existence of molecular bubbles in liquid neon. 
Third, the experimental conditions of MicroCLEAN 
required that the liquid follow the saturation line, render- 
ing it impossible to separate the temperature effects from 
those associated with pressure. It is likely that bubble 
interactions depend on the pressure more strongly than 
on the liquid temperature, and an experiment capable of 
probing the liquid phase space away from the saturation 
line and at higher temperatures could prove interesting. 



VI. CONCLUSION 

In this paper, we have reported measurements of scin- 
tillation light. With PMTs immersed directly in liquid 
neon, we have observed a substantial signal yield of up 



to (3.5 ± 0.4) photoelectrons/keVee. The observed sig- 
nal yield is greater than what would be required for a 
pp-neutrino measurement in CLEAN and much larger 
than the signal yield assumed in previous simulations of 
a tonne-scale CLEAN detector [2,, . We have demon- 
strated the use of activated charcoal as a purifier to re- 
move light reducing impurities in the liquid. We have 
quantified the performance of the prompt ratio PSD 
method for liquid neon by measuring the ERC for a nu- 
clear recoil acceptance probability of approximately 0.5 
between 50 keVr and 300 keVr, although we believe our 
measurement to be an upper limit on the discrimination 
power achievable with liquid neon due to noise in our 
system. As in previous studies of liquid argon, we ob- 
serve a convergence of the fp parameter at low energies. 
We have also made a measurement of CoS in liquid neon. 
Finally, we have observed a very interesting time depen- 
dence of the scintillation of liquid neon, consisting of four 
distinct components and a clear temperature/pressure 
dependence in both intensity and timing. 
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